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Abstract 

We present low-resolution, near-IR JHK spectra of the weak z — 2.39 radio galaxy 53W002, obtained 
with the OH-airglow Suppressor spectrograph (OHS) and Cooled Infrared Spectrograph and Camera for 
OHS (CISCO) on the Subaru Telescope. They cover rest-frame wavelengths of 3400-7200 A, and the 
emission lines of [On] A3727, Rf3, [O in] AA4959, 5007, Ha, [Nn] AA6548, 6583 and [S n] AA6716, 6731 were 
detected. Using the Ha/H/3 line ratio, we find an extinction of E(B — V) — 0.14. The emission- line ratios 
are reproduced by a cloud of electron density n c = 1 x 10 3-4 (cm -3 ) with solar metallicity, ionized by 
an a = —0.7 power-law continuum with ionizing parameter U = 1 X 10~ 3 . In addition to these emission 
lines, we make the first spectroscopic confirmation of the Balmer discontinuity in a high-z radio galaxy. 
Together with rest- frame UV photometry from the literature, we show that at least 1/3 of the present 
stellar mass was formed in the current starburst. The stellar mass was estimated to be (1 — 1.4) x 10 11 Mq 
by one-component model fitting, which is smaller than that of typical z ~ 1 B2/6C radio galaxies. We 
suggest that 53W002 is currently assembling a large part of its stellar mass through merger events with 
the surrounding sub-galactic clumps, some of which can be identified with the La emitters detected in 
narrow-band imaging. After a few such events over the next few Gyr, 53W002 will evolve into a massive 
elliptical galaxy. 
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1. Introduction 

Many high-z field galaxies were discovered in the 1990s 
using the UV drop-out technique. However, these ob- 
jects are rather faint, and infrared spectroscopy is dif- 
ficult (Pettini et al. 1998). On the other hand, high-z 
radio galaxies are relatively bright, and their properties 
have been fairly well-studied. In addition, as the progen- 
itors of present-day massive galaxies, they are important 
for studying galaxy formation. 

53W002 is a z = 2.39 galaxy with relatively weak ra- 
dio power (-51.4GHz ~ 50mJy; Windhorst et al. 1991), 
indicating that its AGN activity is not extreme and it 



1 Based on data collected at Subaru Telescope, which is operated 
by the National Astronomical Observatory of Japan. 



could therefore be suitable for a study of the underlying 
stellar population. Broad-band photometry (Windhorst 
et al. 1991) reveals a possible break at A Tes t ~ 4000 A, 
and suggests recent star-formation activity from spectral 
model fitting. Its rest-frame UV spectrum suggests that 
the main-sequence turnoff age of the current starburst 
is 0.25-0.3 Myr (Windhorst et al. 1991). Radio observa- 
tions have detected lines from CO molecules concentrated 
at the center of the galaxy (Yamada et al. 1995, Scoville 
et al. 1997, Alloin et al. 2000), which may be other indi- 
cations of star-formation activity. 

The radial surface brightness profile of 53W002, as re- 
vealed by HST, is well- modeled by a de Vaucouleurs pro- 
file, suggesting that the galaxy is already dynamically 
relaxed. A V — I radial color gradient suggests that the 
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center of the galaxy consists of a young population (0.3- 
0.5 Gyr), while older stars (0.5-1.0 Gyr) exist in the outer 
(r > 10 kpc) region (Windhorst et al. 1994). 

Previous infrared spectroscopy has revealed the pres- 
ence of the Ha+[Nn] lines in the if-band (Eales, Rawl- 
ings 1993), but a detailed discussion was precluded by 
their limited sensitivity. We have therefore carried out 
low-resolution infrared spectroscopy of 53W002 to inves- 
tigate the properties of its emission-line clouds and un- 
derlying stellar population. For J_ff-band spectroscopy, 
we used a newly-commissioned instrument of the Sub- 
aru telescope, an OH-airglow suppression Spectrograph 
(OHS; Maihara et al. 1994, Iwamuro et al. 2001). OHS 
is a spectroscopic filter which eliminates almost all of 
the OH airglow lines in the J_ff-band, which are the ma- 
jor source of atmospheric background radiation at near- 
infrared wavelengths below 2^im. Using OHS, the back- 
ground level is suppressed by 96% and the sensitivity be- 
comes 1 magnitude fainter. The output from OHS is an 
undispersed beam from which OH airglow is suppressed; 
this is then sent to the Cooled Infrared Camera and Spec- 
trograph for OHS (CISCO; Motohara et al. 1998) to ac- 
quire either images or spectra. CISCO can also be used as 
an ordinary infrared camera and spectrograph attached 
on either the Cassegrain or the Nasmyth focus of the 
telescope. For if-band spectroscopy, we used this mode. 

We describe the observations and data reduction in 
section 2 and present the results in section 3. Discussions 
are given out in section 4 and we summarize our results 
in section 5. We adopt Hq — 65kms _1 and q = 0.1 
throughout this paper. The scale at z — 2.39 is thus 8.5 
kpc arcsec -1 . 

2. Observations and Data Reductions 

2.1. K -Band Spectroscopy 

-ftf-band spectroscopy of 53W002 was carried out on 
1999 May 10, with CISCO mounted at the Cassegrain 
focus of the Subaru telescope. The pixel scale was 
O'.'116/pixcl. A slit width of (f.'7 provided a wavelength 
resolution of ~ 430 at 2.2 fim. The position angle of the 
slit was set to 145° to include object 13 of Pascarelle et 
al. (1998) in the slit. Individual exposures were 200 s, 
and the telescope was nodded slightly (~ 5") after ev- 
ery six frames to provide sky subtraction. A total of 18 
frames were taken for a total exposure time of 3600 s. 
The seeing size was 0'.'4. After the observation, the F0 
star SAO 30082 was observed for correction of the atmo- 
spheric absorption. 

2.2. H-Band Imaging and the JH-Band Spectroscopy 

fZ-band imaging was carried out on 2000 May 20 and 
JiJ-band spectroscopy on 2000 May 20 and 23, with OHS 



and CISCO mounted at the Nasmyth focus of the Subaru 
telescope. The pixel scale was O'.'105/pixel. 

The ii-band images were taken with OHS in the imag- 
ing mode (see Iwamuro et al. 2001). The single expo- 
sure time was 50 s and the telescope was nodded every 3 
frames by 10" in the east-west direction. A total of 12 
frames were acquired, for an exposure time of 600 s. 

The Ji/-band spectra were acquired using a 0'.'95 slit, 
providing a resolution of ~ 210 at 1.65/im. A position 
angle of 90° was used, which was aligned with the major 
axes of the radio and La emission. After each 1000 s 
exposure, the telescope was nodded by 10" along the slit. 
We took 4 frames each day, and the total exposure time 
was 8000 s. The seeing size was 0"6. The F8 star SAO 
30361 and F5 star SAO 30245 were observed after the 
target on May 20 and 23, respectively, for correcting the 
atmospheric absorption. 

2.3. Data Reduction 

The £/-band frames were reduced though standard 
procedures of flat-fielding, sky subtraction, correction of 
bad pixels, residual sky subtraction, before registration 
and coaddition. Photometry was performed in a 1"0 x 2'.'1 
aperture at a position angle of 90°, which is equal to the 
spectroscopic extraction aperture position. The flux in 
the .ff-band was 11.9 ± 0.7 /zJy, calibrated from an obser- 
vation of FS 23. 

The spectroscopic data were reduced in the standard 
manner of flat-ficlding, sky subtraction, correction of bad 
pixels and residual sky subtraction. The spectra of the 
SAO F stars were used to correct for atmospheric extinc- 
tion. 

The wavelength scale of the if-band spectrum was cal- 
ibrated from airglow lines in the raw frames. Such a 
procedure was impossible for the JH spectrum, since all 
the strong OH- lines were suppressed by OHS. We there- 
fore calibrated our spectrum from another observation of 
CISCO without OHS, where OH lines were available. We 
checked the stability of the wavelength scale from all of 
the previous CISCO observations, and found the system- 
atic error to be less than 0.5 pixels (< 5 A). 

Because the seeing size was less than the slit width for 
both spectroscopic observations, we assumed the error in 
the wavelength calibration to be 1/4 of the slit width, 
i.e., ~ ±13 A (1.5 pixels) for the if-band and ~ ±22 A 
(2.4 pixels) for the JH-b&nd. Thus, the systematic error 
in redshifts are 0.0006 in the if -band and 0.0015 in the 
J-ff-band. 

3. Results 

The final JHK spectra are shown in figure 1 . The JH- 
band spectrum was calibrated from the fZ-band photom- 
etry and the if-band spectrum from the K'-band pho- 
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tometry of Yamada et al. (in preparation). A l'.'O x 2'.'1 
aperture, aligned with the position angle used for the 
spectroscopy, was used in both cases. 

3.1. Emission Lines 

We detected emission lines of [O n] A3727, 
H/3, [O in] AA4959, 5007, Ha+[Nn] AA6548, 6583 and 
[Sn] AA6716,6731. [Neiii] A3869 may be present, but this 
region is contaminated with the edge of the Balmer dis- 
continuity. None of the emission lines was resolved at the 
wavelength resolutions of OHS/CISCO, which were 700 
km s _1 in the If -band, and 1500 km s _1 in the /fl-band. 

Due to the low wavelength resolution, the Ha and [N n] 
lines were blended, as was the [Sn] doublet. We de- 
convolved these lines by fitting multiple Gaussian pro- 
files. For Ha+[Nn], we assumed [Nil] A6583/[Nn] A6548 
= 3.0 and demanded that all three lines have the same 
width. For the [Sn] doublet, we fixed the line widths to 
be equal to results from the Ha+[Nn] fitting. The re- 
sults are given in figure 2. The derived line ratios were 
[Nn] A6583/Ha = 2.3±0.6 and [Sn] A6716/[Sn] A6731 = 
0.61 ± 0.44 (2a). Other lines were fitted by single Gaus- 
sian profiles to obtain their redshifts and to determine 
the continuum level; we assumed a common redshift and 
line widths for [O n] , H/3, and [O in] . The emission- line 
properties are listed in table 1 . The fluxes of [O n] , H/3, 
and [O i] are the value above the continuum level. The 
fluxes of the other lines are determined from the results 
of Gaussian fitting, because they are blended. The red- 
shift difference between lines in the JH- and if-bands 
might arise from a systematic error in the wavelength 
calibration. 

The contribution to the broad-band fluxes from the 
emission lines is very large, amounting to 22% in the J- 
band, 30% in the H and 42% in the if '-band. 

We show the spatial intensity profiles of the continuum 
and emission lines along the slit in figure 3. The profiles 
of both [O in] and H«+ [N n] are clearly more extended 
than the adjacent continua. We estimate of the spatial 
extent of the line emission to be roughly 0'.'7-0'.'8, or 6-7 
kpc. 

3.2. Detection of Balmer Discontinuity 

The high sensitivity of OHS allows the continuum of 
53W002 to be clearly detected. We show the SED of the 
continuum in figure 4 with a solid line. We removed all 
of the emission-lines, except for [Nein], using the results 
of the Gaussian profile fitting, and then smoothed the 
resulting spectrum with a 0.05 fxm boxcar filter. 

A clear break in the continuum is seen at 1.3 /zm, which 
is the first spectroscopic detection of the Balmer dis- 
continuity in a high- redshift radio galaxy. The [Nein] 
emission line, which we disregarded above, is located 
here, and may have affected the position of this break. 



To check the effect of this emission line, we sub- 
tract a Gaussian at the position of the line, assuming 
[Nem] A3869/[On] A3727 = 0.2 which was the result of 
simulations from CLOUD Y94 carried out in the following 
analysis. The continuum thus derived is displayed with 
a dotted line in figure 4, and shows a slight shift in the 
position of the break. 

4. Discussion 

4-1. Nature of Emission- Line Clouds 

We measure Ha/H/3 = 3.4. If we adopt an intrinsic 
value of 3.1 (Osterbrock 1989) and the SMC extinction 
curve of Pci (1992) with A v / E(B - V) = 3.1, we derive 
E(B — V) = 0.14. The extinction-corrected emission line 
fluxes are listed in the fifth column of table 1. 

The sulfur line ratio is a good indicator of the electron 
density (Osterbrock 1989), and our measured value of 
[Sn] A6716/[Sn] A6731 = 0.61 ± 0.44 implies n e ~ 1 x 
10 3 cm~ 3 . The 2c lower limit is n c ~ 5 x 10 2 cm -3 , but 
we cannot constrain the upper limit. 

The La flux is 2.8 x 10~ 19 W m" 2 (Keel et al. 1999). 
Because the La morphology is smaller than our slit 
size of l'.'O (Windhorst et al. 1998), we estimate the 
L(La)/L(H/3) ratio to be 3.2. The correction for the 
Galactic extinction of E(B — V)= 0.03 (Burstcin, Hciles 
1982) produces L(La)/L(H/3) = 3.5. A further correc- 
tion for the in situ extinction of E(B — V) = 0.14 results 
in La/H/3 = 17, which is consistent, given the uncer- 
tainty in the extinction, with a low-density value of 23 
(Ferland, Osterbrock 1985). 

In figure 5 we present two emission-line ra- 
tio diagrams. One shows [N n] A6583/Ha versus 
[Om] A5007/H/3 and the other [On] A3727/[Om] A5007 
versus [Om] A5007/H/3. We plot the location of 53W002 
together with the ratios of narrow-line AGNs; H n regions 
taken from the literature. The radio galaxy lies in the re- 
gion occupied by narrow line AGNs and the emission-line 
cloud is therefore ionized by a hidden nucleus. 

To investigate the physical state of the emission-line 
cloud, we carried out photoionization calculations using 
the code CLOUD Y94 (Ferland 1999, 2000). We assume 
a spectral index of a = —0.7 for the ionizing continuum, 
taken from the composite spectra of radio-loud QSOs 
(Cristiani, Vio 1990; Zheng et al. 1997). The observed 
line ratios are reproduced with Z ~ Z Q , \ogU ~ —3, and 
n c = 1 x 10 3 ~ 4 cm~ 3 (see solid and dashed lines in fig- 
ure 5). We could not reproduce the observed ratios with 
an a = —1.5 power law continuum (dot-dashed lines in 
figure 5), nor by models with 0.1 Z & metallicity (dashed- 
long-dashed lines in figure 5). 

We then estimated the mass of ionized gas from the 
physical values derived above and the Ha line luminosity, 
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using the following equations: 

Lua = n 2 c a^ a hv Ha Vf, (1) 
M gas = Vfn e m H , (2) 

where h is the Planck constant, vua the frequency of 
the Ha line, V the volume of the line-emitting cloud, 
/ the filling factor, mn the hydrogen mass, and a^ a — 
6.04 x 10- 14 (cm 3 s- 1 ) (Osterbrock 1989) the Ha recom- 
bination coefficient under Case B. We derive a gas mass 
of 10 7 - 8 M o . 

A rough estimate of the filling factor, /, of the cloud 
can be made by assuming the cloud to be a sphere of 3 
kpc radius (section 3.1). We calculate / ~ 10~ 5 -10~ 7 , 
consistent with the typical values of low- z radio galaxies 
of 10- 4 -10 -6 (Heckman et al. 1982; van Breugel et al. 
1985). 

High-resolution HST imaging shows the presence of an 
extended blue cloud which peaks 0'/45 west of the nucleus 
(Windhorst et al. 1998). Our emission- line region, how- 
ever, almost coincides with the continuum core. There- 
fore, what we are seeing is not those blue clouds, but a 
cloud which envelops the whole galaxy. This gas could 
be a remnant of a previous merger event, resulting in the 
starburst activity observed in our spectra as the Balmcr 
discontinuity (see section 4.2.2). Metallicity as high as 
Z Q might be produced during this period. 

4-2. SED and Age of the Host Galaxy 

The spectral energy distribution of 53W002 is shown 
by the open circles and squares in figure 6. The optical 
data were taken from Windhorst et al. (1991) and the 
infrared data from line-free regions of our spectra. An 
additional data point was derived from the continuum 
after subtracting the [Ne in] line assuming [Ne m]/[0 n] 
= 0.2±0.1. Because Windhorst et al. (1991) measured 
the total magnitudes, our infrared data were multiplied 
by a factor of 1.2 to match their H- and K-h&nd pho- 
tometry. We disregard the U -band data of Windhorst et 
al. (1991), which is shortward of La in the rest-frame, 
and is likely to be strongly attenuated by the La forest 
from the surrounding objects. 

Next, we removed two non-stellar components from the 
SED. One is the nebular continuum emission and the 
other is the scattered light from the hidden nucleus. 

The strength of the nebular continuum was calculated 
from the H/3 emission line flux using the tabulated co- 
efficients of Abler (1984) for T c = 1 x 10 4 K. We ig- 
nore the contribution from He recombination, which our 
CLOUDY94 simulations indicate will be weak. 

We assume the intrinsic shape of the continuum from 
the hidden nucleus to be an a = —0.7 power law. Since 
wc clo not know whether the scatterers are dust grains 
or electrons, we consider both cases. Because electron- 
scattering is achromatic, the scattered spectrum is the 
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same as the incident spectrum, while the dust-scattered 
spectrum was from Cimatti et al. (1994). After applying 
the extinction of E(B — V) = 0.14, the flux was scaled 
to match the upper limits of the scattered AGN contri- 
bution estimated by Windhorst et al. (1998) as shown in 
figure 6. 

The resulting SEDs of the underlying stellar popula- 
tion, after subtracting these two components, is shown 
by the filled circles and triangles in figure 6. We tried to 
fit these SEDs with a few models of galactic spectra. 

4-2.1. One- component model 

First, we fitted these SEDs with single-component 
synthetic spectra produced by the spectrophotomet- 
ric galaxy evolution code PEGASE (Fioc, Rocca- 
Volmcrange 1997), to which an extinction of E(B — V) = 
0.14 was applied. Instead of the SMC extinction curve, 
we applied the extinction curve of Calzetti et al. (1994), 
which is an empirical relation for the stellar UV con- 
tinuum derived from the spectra of nearby starburst re- 
gions. We considered three star-formation histories - 
an instantaneous burst, and two exponential bursts with 
e-folding times of 200 Myr and 1 Gyr — and the age was 
left as a free parameter. We adopted the Scalo IMF. We 
also fit our models to the observed SED, to investigate 
the effect of subtracting the non-stellar emission. These 
results are given in table 2. 

The results are not very sensitive to the assumed non- 
stellar contributions, since the age of the galaxy is mainly 
determined by the strength of the Balmcr discontinu- 
ity. Although the fitted parameters vary according to 
the assumed star-formation history, the age is always 
young: 70-700 Myr, corresponding to a formation red- 
shift Zf orm = 2.5-3.2. 

Windhorst et al. (1991) showed, from the far-UV SED 
and possible detections of UV absorption lines, that the 
main-sequence cutoff age is 250-320 Myr, so we adopt 
the r = 200 Myr model in the following discussions. The 
stellar mass of the galaxy is then ~ 10 11 Mq, which is of 
the order of that of an L* galaxy. 

4-2.2. Two- component model 

According to the high-resolution HST imaging, 
53W002 has a radial color gradient, suggesting that the 
outskirts of the galaxy is comprised of older stars (0.5-1 
Gyr; Windhorst et al. 1994). There may therefore have 
been additional star-formation episodes prior to the cur- 
rent one, and we attempt to fit the SED with two stellar 
components to evaluate how much stellar mass was cre- 
ated by the current starburst activity. 

The models consisted of an old component of fixed age 
and a younger component of unknown age. A t = 200 
Myr exponential burst with Scalo IMF was assumed for 
both components. We varied the fraction of the initial 
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gas mass in the old component from 0.1 to 0.9 in 0.1 steps. 
Two ages for the old component were assumed: one is 1 
Gyr, the age suggested by Windhorst et al. (1994); the 
other is 3 Gyr, the maximum age allowed by the age of 
the universe. The results of the model fitting arc shown 
in table 3 and figure 7. 

The x 2- test was used to assess the results of the fit- 
ting. We reject models which are not consistent with the 
data at the 99% confidence level, which corresponds to 
X 2 > 25. We therefore exclude models where the mass 
fraction of old component is 0.9 or larger. We also reject 
models which produce an age of less than 200 Myr for 
the current burst, which is a lower limit derived from the 
UV spectrum (Windhorst et al. 1991), and only those 
models whose old mass fraction is less than 0.7 remain. 
We therefore conclude that 53W002 is undergoing a mas- 
sive starburst, and that at least 1/3 of its stellar mass is 
forming in the current burst. 

4.2.3. Mass and formation of 53W002 

We next compare the stellar mass with those of other 
radio galaxies in the literatures. Because these mass es- 
timations were carried out by single-component SED fit- 
ting, we use the mass of 53W002 derived using the one- 
component model. 

The stellar mass of the galaxy was estimated to be less 
than 1.4 x 10 11 M by the one-component model, similar 
to that of z = 1.824 radio galaxy 3C 256 (Simpson et 
al. 1999). This galaxy also underwent a recent starburst 
and its stellar mass is only 1.2 x 10 11 M© (converted to 
our adopted cosmology). Best et al. (1998) estimated the 
masses of 26 z ~ 1 3C radio galaxies using the GISSEL 
spectral synthesis codes of Bruzual and Chariot (1993), 
and derived a median value of 4.4± 1.7 x 10 11 M (again, 
converted to our adopted cosmology) . We determine that 
the effect of the different IMFs used is only 10%, and 
therefore conclude that the mass of 53W002 is less than 
1/3 of a typical z ~ 1 3C radio galaxy. However, Eales 
et al. (1997) reported a correlation between the radio 
and i-T-band luminosities of z ~ 1 radio galaxies, and 
the radio luminosity of 53W002 is similar to those of the 
B2/6C radio galaxies, which are on average ~ 0.6 mag 
fainter than the 3C objects. If we assume that the stellar 
mass correlates with the X-band luminosity, the mass of 
z ~ 1 B2/6C radio galaxies is ~ 2.5 x 10 11 M , which is 
still twice the value of 53W002. 

There are 17 candidate La-emitters in the field of 
53W002, and 7 of them have been spectroscopically con- 
firmed to lie at the same redshift as the radio galaxy 
(Pascarelle et al. 1996a,b,1998; Keel et al. 1999). These 
objects are of sub-galactic size, and were postulated to be 
'building blocks' which will merge into a luminous galaxy 
by the present epoch. If 53W002 experiences additional 
starburst activity every 500 Myr by merging with these 



surrounding sub-galactic clumps, it will grow into a mas- 
sive galaxy of a fewxlO n M0 by z = 1. These results 
strongly support the hypothesis that 53W002 is a form- 
ing massive elliptical galaxy undergoing multiple merger 
events, and an ancestor of present-day luminous galaxies. 

5. Summary 

We have performed infrared JHK spectroscopy of the 
weak radio galaxy 53W002 at z = 2.39, using CISCO and 
the newly commissioned instrument OHS on the Subaru 
telescope. 

We detected the rest-frame optical emission lines of 
[0 11] A3727, H/3, [O m] AA4959, 5007, 

Ha, [Nil] AA6548, 6583 and [S 11] AA6716, 6731. These 
emission lines were reproduced by a model cloud with 
Z ~ Z Q , \ogU ~ -3, and n c = 1 x 10 3 - 4 (cuT 3 ), 
ionized by an a = —0.7 power-law continuum. The 
mass of ionized gas in the cloud was estimated to be 
-Mgas = 10 7 ~ 8 M Q , and the filling factor was calculated to 
be as small as 10~ 5 -10~ 7 , indicating significant dumpi- 
ness. We found no offset between the peak positions of 
the line and continuum emission, which suggests that the 
emission line cloud is not the extended blue cloud re- 
ported in Windhorst et al. (1998). 

The Balmer discontinuity was detected clearly for the 
first time in the continuum of a high-rcdshift radio 
galaxy. We were able to fit the rest-frame UV-optical 
continuum by a 70-700 Myr-old stellar population. A fit 
with two-component model indicated that at least one- 
third of the total stellar mass is involved in the current 
starburst. We estimated the total stellar mass to be 1- 
1.4 x 10 11 Mq by one-component model, which is smaller 
than typical z ~ 1 B2/6C radio galaxies. The mass be- 
comes twice as much if we introduce the two-component 
model, but is still smaller than that of z ~ 1 3C radio 
galaxies. We therefore conclude that 53W002 is now as- 
sembling a large part of its stellar mass through mergers 
with the surrounding sub-galactic objects, and will end 
up as a massive elliptical observed today. 

We are indebted to all staff members of the Subaru 
telescope, NAOJ. We appreciate M. Fioc and B. Rocca- 
Volmerange for generously offering their galaxy modeling 
code, PEGASE and G. Ferland for the spectral synthesis 
code CLOUDY94. 
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Fig. 1. JHK spectra of 53W002. The lower trace shows the lcr noise level calculated from the background level. Many 
strong forbidden lines and weak H/3 line are present. Due to the low-wavelength resolution, Ha+[Nll] AA6548,6583 and 
[S n] AA6716,6731 are blended. 




Fig. 2. Deconvolution of the blended Ha+[Nll] and [S II] lines. 
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Fig. 3. Spatial profiles of the spectra of 53W002. In the upper panels, west is to the right, while south-east is to the right 
in the lower two. The dotted lines arc fitted Gaussian profiles. Both [O in] and Ha+[N II] lines are more spatially 
extended than the continuum. The seeing sizes for the //-band and iiT-band spectra are 0'.'6 and 0'.'4, respectively. 



10 20 30 40 50 10 20 30 40 50 




10 20 30 40 10 20 30 40 

Pixels 



Fig. 4. Continuum of 53W002. All of the emission lines, except [Ncm], have been subtracted by Gaussian fitting, and the 
resulting spectrum has been smoothed with a 0.05 fim boxcar filter. The solid line shows the raw continuum, while the 
dotted line shows further subtraction of [Ncm], where we have assumed [Ne m]/[0 II] = 0.2 (see text). 
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Fig. 5. Reddening corrected line ratios of 53W002 (open circle with error bars) plotted together with those of Seyfert 
2s (filled circles), radio galaxies (filled squares), Hll regions (open circles), Hn galaxies (open circles), and starburst 
galaxies (open triangles) taken from the literature (Veilleux, Ostcrbrock 1987; McCall et al. 1985; French 1980; Koski 
1978; Costero, Ostcrbrock 1977). The results of photoionization calculations carried out using CLOUDY94 arc indicated 
by the lines. Solar metallicity and a spectral index of a = —0.7 for the ionizing continuum were assumed, except for 
one case where we assumed a = —1.5 (thin dot-dashed line), and one case of Z = O.IZq (thin dash-long-dashcd line). 
The ionization parameter (U) varies along each line, with representative points labeled by the value of log U. 




Fig. 6. Spectral energy distributions of 53W002. Open circles and triangles are observed values from the present observations 
and Windhorst et al. (1991), respectively. The assumed AGN and nebular continuum components were subtracted (see 
text) to produce the SED of the stellar component alone, which is shown by the filled symbols. The filled squares 
indicate the upper limit of the unresolved AGN component determined from HST observations (Windhorst et al. 1998). 
The nebular continuum was calculated from the H/3 flux, and is shown by the long dashed line. The solid lines show 
the best-fit one-component spectral model. The dashed line shows the AGN contribution, which is assumed to be 
dust-scattered in the left panel, and electron-scattered in the right panel. 
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Fig. 7. Results of fitting two-component models to the stellar SED of 53W002. The AGN component is assumed to be a 
dust-scattered power-law. The dotted lines show the spectra of the old component, while the dashed lines show that of 
the current burst. 
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Table 1. Emission line properties of 53W002. All errors are 2a. 



Line 


Arest 


Redshift* Flux 


Flux (corrected) 
(xl0~ 19 Wm~ 2 ) 


Continuum 


EW last 




(A) 


(xl0~ 19 Wm~ 2 ) 


(xlO" 18 Wm-Vm) 


(A) 


[On] 


3727 


2.3923±0.0004 4.8±0.4 


8.8±0.8 


6.2±1.3 


230±50 


H/3 


4861 


2.3923±0.0004 0.88±0.18 


1.4±0.3 


10.3±0.8 


20±8 


[Ora] 


4959 


2.3923±0.0004 2.1±0.3 


3.3±0.4 


10.3±0.9 


61±9 


[Om] 


5007 


2.3923±0.0004 7.3±0.5 


11.3±0.8 


10.3±0.9 


210±20 


[Oi] 


6300 


2.3833±0.0076 0.76±0.36 


1.1±0.5 


5.9±1.4 


37±20 


Ha 


6563 


2.3895±0.0008 3.1±0.7 


4.3±0.9 


5.0±1.6 


180±70 


[Nil] 


6583 


2.3895±0.0008 7.2±0.8 


9.9±1.1 


5.0±1.6 


420±140 


[S ii] 


6716 


2.3890±0.0032 1.0±0.7 


1.4±1.0 


5.2±2.5 


60±60 


[Sn] 


6731 


2.3890±0.0032 1.6±0.7 


2.2±1.0 


5.2±2.5 


90±70 



* No systematic errors are included. 



Table 2. Parameters for the best-fit models of the spectral energy distribution. 



Galaxy evolution model Non-stellar* x 2 Age (Myr) Stellar mass (M Q ) SFR (M© yr L ) 



Instantaneous burst . . . 




N 


19.1 


80 


9.7 x 


10 iu 









E 


19.0 


70 


6.8 x 


10 10 









D 


19.1 


80 


8.1 x 


10 10 





Exponential burst (r = 


200 Myr) 


N 


21.3 


350 


1.3 x 


10 11 


1.5 x 10 2 






E 


12.5 


350 


9.7 x 


10 10 


1.1 x 10 2 






D 


10.0 


350 


1.1 X 


10 11 


1.2 x 10 2 


Exponential burst (r = 


1 Gyr) 


N 


23.8 


600 


1.4 x 


10 11 


1.8 x 10 2 






E 


11.0 


600 


1.1 X 


10 11 


1.4 x 10 2 






D 


8.9 


700 


1.3 x 


10 11 


1.4 x 10 2 



* Model of non-stellar components. N: No non-stellar component is assumed. D: Nebular continuum + dust scattered 
a = —0.7 power law. E: Nebular continuum + electron scattered a — —0.7 power law. 
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Table 3. Best-fit two-components models of the spectral energy distribution. 



Age of old component 




1 Gy 


r 




3 Gy 


r 


Non-stellar* Old frac. T 


x 2 


Age* 
(Myr) 


Stellar mass 
(M e ) 


x 2 


Age* 
(Myr) 


Stellar mass 
(M ) 


D 0.1 


9.5 


350 


1.5x10" 


9.2 


350 


1.2x10" 


0.2 


11.4 


350 


1.6x10" 


9.0 


350 


1.3x10" 


0.3 


8.9 


300 


1.6xl0 n 


9.9 


350 


1.5x10" 


0.4 


10.0 


300 


1.7x10" 


9.6 


300 


1.5x10" 


0.5 


8.6 


250 


1.8x10" 


8.4 


300 


1.7x10" 


0.6 


9.5 


200 


1.9x10" 


11.8 


300 


2.1x10" 


0.7 


9.8 


160 


2.1x10" 


12.2 


250 


2.3x10" 


0.8 


13.3 


80 


2.4x10" 


18.7 


180 


2.7x10" 


0.9 


84.4 


40 


2.9x10" 


45.8 


70 


3.7x10" 


E 0.1 


12.4 


300 


1.3x10" 


13.5 


300 


9.3xl0 lu 


0.2 


11.2 


300 


1.4x10" 


12.1 


300 


1.0x10" 


0.3 


12.8 


300 


1.5x10" 


11.3 


300 


1.2x10" 


0.4 


11.7 


250 


1.5x10" 


11.7 


300 


1.4x10" 


0.5 


14.2 


250 


1.6x10" 


14.6 


300 


1.6x10" 


0.6 


13.1 


180 


1.7x10" 


13.8 


250 


1.7x10" 


0.7 


14.6 


140 


1.9x10" 


18.1 


200 


1.9x10" 


0.8 


19.5 


50 


2.2x10" 


24.4 


160 


2.4x10" 


0.9 


103.4 


40 


2.6x10" 


54.5 


60 


3.4x10" 



* Model of non-stellar components. D: Nebular continuum + dust scattered a = —0.7 power law. E: Nebular continuum + 
electron scattered a = —0.7 power law. 
f Mass fraction of old component. 
J Age of current burst component. 



